Copy  ji£,  of  ST)  C'opie 


U.  S.  Navy  Underwater  Sound  Laboratory 
Fort  Trumbull,  New  London,  Connecticut 


INTENSITY  J^CULATIONS  J^pING  A 
GEOMETRICAL  ^ODEL  OF  SifEIVS  LAW 


PISTRIBUT2 * * * * 7  ON  STATEMENT l 

Approved  for  public  release; 
Distribution  Unlimited 


HTBODUCTION 


(l)  The  formula  for  Io/l  here  used  does  not  take  into  account  any 
intensity  loss  that  may  be  incurred  by  attenuation.  Spreading  loss  alone 
is  considered. 


(2)  35ie  spreading  loss  formula  obtained  here  is  invalid  for  calcula- 

tions of  intensity  at  all  points  on  the  limit  eg  ray  beyond  the  point  where 

the  ray  reaches  its  vertex  velocity. 


TERMINOLOGY  AND  BASIC  EQ1 ATIONS 


Before  going  any  further  the  matter  of  terminology  haB  to  be  somewhat 

better  organized.  In  general  the  velocity  of  sound  varies  continuously 
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with  depth  and  this  variation  may  not  be  llnea-.  In  refraction  calcula- 
tions It  Is  therefore  necessary  to  approximate  the  nonlinear  variation  of 
velocity  with  something  linear,  thus  creating  ayero  of  fixed  velocity 
gradients.  In  the  discussion  that  follows,  those  layers  will  be  Indicated 
by  the  symbol  n (n  » 1*2, 3, 4, . . . ) and  numbered  consecutively  along  the 
ray  path  in  such  a manner  that  the  number  of  the  layer  increases  by  one 
each  time  the  ray  crosses  a layer  boundary  or  is  reflected  from  it.  For 
the  purposes  of  simplified  expression  the  point  of  origin  of  each  ray,  as 
well  as  the  point  of  its  termination,  (the  latter  a matter  of  choice  in 
each  calculation)  are  considered  as  layer  boun-laries  (see  Fig.  l). 


Figure  1 

When  n is  attached  as  a subscript  to  stJie  other  symbol,  it  indicates 
that  particular  value  of  the  symbol  appropriate  at  the  point  where  the 
sound  ray  crosses  the  boundary  between  layer  (n-l't  and  layer  n (see 

Fig.  1). 

By  way  of  completing  our  picture  of  the  \ elocity  structure  of  the 
ocean  it  is  convenient  to  associate  with  each  layer  a constant  velocity 

gradient: 
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C sec-*) 


where: 


V^S-  the  velocity  of  Bound  at  the  points  indicated  (ft/eec) 
^n*S-  the  depth  at  the  same  points  (ft) 


As  a starting  point  the  following  equation  will  be  used  (for  deri- 
vation see  "Fundamentals  of  Sonar,"  J.  W.  Horten,  U.  S.  Naval  Institute, 
1957,  pages  99-100): 


where:  0's  - indicate  the  angles  of  inclination  (with  respect  to  the 

horizontal)  of  the  ray  (radians  or  degrees) 

N - total  number  of  layers  passed  through  by  ray 

R - total  horizontal  range  (yds) 

B»e  absolute  value  signs  are  utilized  because  w*  are  interested  in  the 
separation  between  the  bounding  rays  of  a ray  bundle  and  not  in  the 
direction  in  which  this  separation  is  measured.  This  equation  is  valid 
no  matter  what  path  the  ray  follows  between  its  origin  and  the  terminal 
points.  To  express  it  in  terms  of  the  geometric  model  we  also  need  the 
relations  between  the  trigonometric  functions  involved  and  the  geometric 
model.  The  basic  relation  here  is  Snell 'a  Lav: 

(3)  - -siss;  (i^) 


where  Vx  is  the  vertexing  velocity  of  the  sow  d ray.  As  the  Vn's 
ar*  fixed  by  the  velocity  structure  of  the  ocean,  for  a fixed  source 
the  only  Independent  variable  is  ©i  . For  eacl  a different  ray 

path  is  specified  because  the  ratio  Vx  is  fixed  by  Snell’s  Law,  thus 
making  all  ©n's  (n^l)  dependent  on  the  choice  cf  . Assuming  thus  a 
f^*®^  path  with  known  Vx  we  have  the  following  two  relations: 
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Slrt«„  ~ * (+  vhen  Oq  is  positive;  - when  On  Is 

negative).  Thin  rule  is  reversed  at 
each  reflection  from  what  it  was  before 
that  reflection. 

These  relations  vill  enable  us  to  obtain  the  f.nal  expression  for  Iq/I 
entirely  in  terms  of  the  velocity  structure  for  each  ray  of  fixed  V*  . 

Before  this  we  have  two  ether  factors  (R  tnd  ) involving  the  total 
range  R which  must  be  determined.  In  this  dr-termination  we  «han  use 
two  generalizations  from  the  incremental  range  equations  presented  in  USL 
Memorandum  serial  1110-83 . They  are: 


where : Rq's  - Indicate  the  horizontal  Incremental  ranges  covered  by 
sound  ray  in  each  constant  gradient  layer  (yds). 

Ifce  second  form  of  the  above  equation  (5b)  is  tie  result  of  direct  sub- 
stitution from  (4b)  and  the  sign  convention  of  the  latter  is  retained.  It 
is  important  to  remember  that  in  ray  calculations  a ray  directed  downward 
forms  a positive  angle  with  the  horizontal  and  the  ray  directed  upwards  a 
negative  angle. 

The  derivative  dR/d0i  is  calculated  as  shown  in  the  appendix  using 
(5»)  sad  (3).  Then,  by  means  of  (4a)  and  (4b)  She  derivative  is  expressed 
in  terms  of  the  velocity  structure. 


USL  Tech.  Memo. 

No.  1110-44-58 


GENERAL  CASE  (N  LAYERS) 


,,  „ general  expression  is  obtained  now  by  "ibstitution  from  Eqs.  (4a). 
('♦b),  (5a)  and  (A7)  into  Eq.  (2)  ^ " 

L = vjiLRn  \lfa R* 

j v,  vr  Z.< 


Since  the  velocities  are  positive  and  the  two  s iuare  roots  may  be  always 
regarded  as  positive  we  can  write  the  final  fora  as  follows: 


(«>  -i»  = 


v;  -v,‘  /v.’  -vi 


!>•) 


(tVyF^v* 


As  stated  before  this  equation  is  applieab  .e  in  all  cases  except  the 
tW°  ^Sinning.  in  all  cases  conf  ision  will  be  avoided  if 

Eq.  1,6;  is  reduced  to  its  simplest  form  for  eac  i specific  case  before  any 
numerical  calculations  are  made. 

A special  case,  sometimes  giving  rise  to  c >nfusion,  concerns  the 
validity  of  the  equation  after  a ray  has  passed  its  vertex ing  point  and 
is  crossing  the  layers  in  the  opposite  directioi  to  that  preceding  the 
vertexlng  point.  In  this  case  the  sign  rules  should  be  used  as  usual, 
taking  into  consideration  any  reflection?  that  any  have  occurred,  where 
they  apply  in  Eq.  (6).  The  choice  of  lay  ers  an  l layer  boundaries  should 
also  be  made  in  accordance  with  the  directions  previously  given.  Spe- 
cifically, this  means  that  the  point  of  vertex!  ig  is  not  considered  as  a 

“f  thQt  thC  velocity  Vn  r Vnf,  (n  designating  vertexing 
layer;  but  the  exit  and  entrance  angles  are  opposite  in  sign.  The  range 
ncrement  too  in  this  case  would  be  the  distance  between  the  two  points 
at  which  the  ray  crosses  the  same  boundary  in  opposite  directions  (see 
Fig.  1). 


0.  RJRA'S 
Physic!  lit 
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In  Eq.  (Al)both  Qa  and  ©n  ^ j are  still  dependent  on  O3.  and  this  must 
be  taken  account  of  when  differentiating.  Th.'  s operation  is  simplified, 
however,  by  avoiding  the  immediate  replacement  of  these  terms  by  functions 
Involving  Cl  . Performing  the  differentiation  we  obtain: 
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